Collaborative Expression: Transcriptomics of Conus virgo Suggests Contribution of Multiple Secretory Glands to Venom Production.

Venomous marine gastropods of the family Conidae are among the most diversified predators in marine realm-in large due to their complex venoms. Besides being a valuable source of bioactive neuropeptides conotoxins, cone-snails venoms are an excellent model for molecular evolution studies, addressing origin of key innovations. However, these studies are handicapped by scarce current knowledge on the tissues involved in venom production, as it is generally assumed the sole prerogative of the venom gland (VG). The role of other secretory glands that are present in all Conus species (salivary gland, SG) or only in some species (accessory salivary gland, ASG) remains poorly understood. Here, for the first time, we carry out a detailed analysis of the VG, SG, and ASG transcriptomes in the vermivorous Conus virgo. We detect multiple transcripts clusters in both the SG and ASG, whose annotations imply venom-related functions. Despite the subsets of transcripts highly-expressed in the VG, SG, and ASG being very distinct, SG expresses an L-, and ASG-Cerm08-, and MEFRR- superfamily conotoxins, all previously considered specific for VG. We corroborate our results with the analysis of published SG and VG transcriptomes from unrelated fish-hunting C. geographus, and C. striatus, possibly fish-hunting C. rolani, and worm-hunting Conus quercinus. In spite of low expression levels of conotoxins, some other specific clusters of putative venom-related peptides are present and may be highly expressed in the SG of these species. Further functional studies are necessary to determine the role that these peptides play in envenomation. In the meantime, our results show importance of routine multi-tissue sampling both for accurate interpretation of tissue-specific venom composition in cone-snails, and for better understanding origin and evolution of venom peptides genes.

Oxidative Folding Catalysts of Conotoxins Derived from the Venom Duct Transcriptome of C. frigidus and C. amadis.

Two novel redox conopeptides with proline residues outside and within the active site disulfide loop were derived from the venom duct transcriptome of the marine cone snails Conus frigidus and Conus amadis. Mature peptides with possible post-translational modification of 4-trans-hydroxylation of proline, namely, Fr874, Fr890[P1O], Fr890[P2O], Fr906, Am1038, and Am1054, have been chemically synthesized and characterized using mass spectrometry. The estimated reduction potential of cysteine disulfides of synthetic peptides varied from -298 to -328 mV, similar to the active site cysteine disulfide motifs of the redox family of proteins. Fr906/Am1054 exhibited pronounced catalytic activity and assisted in improving the yields of natively folded globular form α-conotoxin ImI. Three-dimensional (3D) structures of the redox conopeptides were optimized using computational methods and verified by 2D-ROESY NMR spectroscopy: C. frigidus peptides adopt an N-terminal helical fold and C. amadis peptides adopt distinct structures based on the Phe4-Pro/Hyp5 peptide bond configuration. The shift in the cis-trans configuration of the Phe4-Pro/Hyp5 peptide bond of Am1038/Am1054 was observed between reduced free thiol and oxidized disulfide forms of the optimized peptides. The report confirms the position-specific effect of hydroxyproline on the oxidative folding of conotoxins and sequence diversity of redox conopeptides in the venom duct of cone snails.

Systematic dissection of genomic features determining the vast diversity of conotoxins.

BACKGROUND: Conus, a highly diverse species of venomous predators, has attracted significant attention in neuroscience and new drug development due to their rich collection of neuroactive peptides called conotoxins. Recent advancements in transcriptome, proteome, and genome analyses have facilitated the identification of conotoxins within Conus' venom glands, providing insights into the genetic features and evolutionary patterns of conotoxin genes. However, the underlying mechanism behind the extraordinary hypervariability of conotoxins remains largely unknown. RESULTS: We analyzed the transcriptomes of 34 Conus species, examining various tissues such as the venom duct, venom bulb, and salivary gland, leading to the identification of conotoxin genes. Genetic variation analysis revealed that a subset of these genes (15.78% of the total) in Conus species underwent positive selection (Ka/Ks > 1, p < 0.01). Additionally, we reassembled and annotated the genome of C. betulinus, uncovering 221 conotoxin-encoding genes. These genes primarily consisted of three exons, with a significant portion showing high transcriptional activity in the venom ducts. Importantly, the flanking regions and adjacent introns of conotoxin genes exhibited a higher prevalence of transposon elements, suggesting their potential contribution to the extensive variability observed in conotoxins. Furthermore, we detected genome duplication in C. betulinus, which likely contributed to the expansion of conotoxin gene numbers. Interestingly, our study also provided evidence of introgression among Conus species, indicating that interspecies hybridization may have played a role in shaping the evolution of diverse conotoxin genes. CONCLUSIONS: This study highlights the impact of adaptive evolution and introgressive hybridization on the genetic diversity of conotoxin genes and the evolution of Conus. We also propose a hypothesis suggesting that transposable elements might significantly contribute to the remarkable diversity observed in conotoxins. These findings not only enhance our understanding of peptide genetic diversity but also present a novel approach for peptide bioengineering.

The role of the Sunda shelf biogeographic barrier in the cryptic differentiation of Conus litteratus (Gastropoda: Conidae) across the Indo-Pacific region.

Geographical and oceanographic processes have influenced the speciation of marine organisms. Cone snails are marine mollusks that show high levels of endemism and a wide distributional range across the Indian and Pacific Oceans. Discontinuities in distributions caused by biogeographic barriers can affect genetic connectivity. Here we analysed the connectivity within Conus litteratus using samples from the Lakshadweep archipelago (Arabian Sea, Indian Ocean) and from the Pacific Ocean. Maximum likelihood analyses based on the mitochondrial cytochrome C oxidase subunit I (COI) and on the non-coding 16S ribosomal RNA (16S rRNA) genes revealed cryptic diversity within C. literatus occupying distinct oceanographic regions. The intraspecific genetic distances between the two distinct clades of C. literatus from the Arabian Sea and the Pacific Ocean ranged from 7.4% to 7.6% for COI and from 2.4% to 2.8% for 16S rRNA genes, which is larger than the threshold limit for interspecific differentiation. The haplotype network analysis also corroborated the existence of two different lineages within C. litteratus. The detected genetic discontinuities reflect the effect of the Sunda shelf biogeographic barrier on the allopatric divergence of C. litteratus.

Mitogenome Characterization of Four Conus Species and Comparative Analysis.

Cone snails, as a type of marine organism, have rich species diversity. Traditionally, classifications of cone snails were based mostly on radula, shell, and anatomical characters. Because of these phenotypic features' high population variability and propensity for local adaptation and convergence, identifying species can be difficult and occasionally inaccurate. In addition, mitochondrial genomes contain high phylogenetic information, so complete mitogenomes have been increasingly employed for inferring molecular phylogeny. To enrich the mitogenomic database of cone snails (Caenogastropoda: Conidae), mitogenomes of four Conus species, i.e., C. imperialis (15,505 bp), C. literatus (15,569 bp), C. virgo (15,594 bp), and C. marmoreus (15,579 bp), were characterized and compared. All 4 of these mitogenomes included 13 protein-coding genes, 2 ribosomal RNA genes, 22 tRNA genes, and non-coding regions. All the Protein Codon Genes (PCGs) of both newly sequenced mitogenomes used TAA or TAG as a terminal codon. Most PCGs used conventional start codon ATG, but an alternative initiation codon GTG was detected in a gene (NADH dehydrogenase subunit 4 (nad4)) of C. imperialis. In addition, the phylogenetic relationships were reconstructed among 20 Conus species on the basis of PCGs, COX1, and the complete mitogenome using both Bayesian Inference (BI) and Maximum Likelihood (ML). The phylogenetic results supported that C. litteratus, C. quercinus, and C. virgo were clustered together as a sister group (PP = 1, BS = 99), but they did not support the phylogenetic relation of C. imperialis and C. tribblei (PP = 0.79, BS = 50). In addition, our study established that PCGs and complete mitogenome are the two useful markers for phylogenetic inference of Conus species. These results enriched the data of the cone snail's mitochondrion in the South China Sea and provided a reliable basis for the interpretation of the phylogenetic relationship of the cone snail based on the mitochondrial genome.

Anti-Ovarian Cancer Conotoxins Identified from Conus Venom.

Conotoxins constitute a treasury of drug resources and have attracted widespread attention. In order to explore biological candidates from the marine cone snail, we isolated and identified three novel conopeptides named as Vi14b, Vi002, Vi003, three conotoxin variants named as Mr3d.1, Mr3e.1, Tx3a.1, and three known conotoxins (Vi15a, Mr3.8 and TCP) from crude venoms of Conus virgo, Conus marmoreus and Conus texile. Mr3.8 (I-V, II-VI, III-IV) and Tx3a.1 (I-III, II-VI, IV-V) both showed a novel pattern of disulfide connectivity, different from that previously established for the µ- and ψ-conotoxins. Concerning the effect on voltage-gated sodium channels, Mr3e.1, Mr3.8, Tx3a.1, TCP inhibited Nav1.4 or Nav1.8 by 21.51~24.32% of currents at semi-activated state (TP2) at 10 µmol/L. Certain anti-ovarian cancer effects on ID-8 cells were exhibited by Tx3a.1, Mr3e.1 and Vi14b with IC50 values of 24.29 µM, 54.97 µM and 111.6 µM, respectively. This work highlights the role of conotoxin libraries in subsequent drug discovery for ovarian cancer treatment.

Vexitoxins: conotoxin-like venom peptides from predatory gastropods of the genus Vexillum.

Venoms of predatory marine cone snails are intensely studied because of the biomedical applications of the neuropeptides that they contain, termed conotoxins. Meanwhile some gastropod lineages have independently acquired secretory glands strikingly similar to the venom gland of cone snails, suggesting that they possess similar venoms. Here we focus on the most diversified of these clades, the genus Vexillum. Based on the analysis of a multi-species proteo-transcriptomic dataset, we show that Vexillum species indeed produce complex venoms dominated by highly diversified short cysteine-rich peptides, vexitoxins. Vexitoxins possess the same precursor organization, display overlapping cysteine frameworks and share several common post-translational modifications with conotoxins. Some vexitoxins show sequence similarity to conotoxins and adopt similar domain conformations, including a pharmacologically relevant inhibitory cysteine knot motif. The Vexillum envenomation gland (gL) is a notably more recent evolutionary novelty than the conoidean venom gland. Thus, we hypothesize lower divergence between vexitoxin genes, and their ancestral 'somatic' counterparts compared to that in conotoxins, and we find support for this hypothesis in the evolution of the vexitoxin cluster V027. We use this example to discuss how future studies on vexitoxins can inform the origin of conotoxins, and how they may help to address outstanding questions in venom evolution.

Reconstructing the Origins of the Somatostatin and Allatostatin-C Signaling Systems Using the Accelerated Evolution of Biodiverse Cone Snail Toxins.

Somatostatin and its related peptides (SSRPs) form an important family of hormones with diverse physiological roles. The ubiquitous presence of SSRPs in vertebrates and several invertebrate deuterostomes suggests an ancient origin of the SSRP signaling system. However, the existence of SSRP genes outside of deuterostomes has not been established, and the evolutionary history of this signaling system remains poorly understood. Our recent discovery of SSRP-like toxins (consomatins) in venomous marine cone snails (Conus) suggested the presence of a related signaling system in mollusks and potentially other protostomes. Here, we identify the molluscan SSRP-like signaling gene that gave rise to the consomatin family. Following recruitment into venom, consomatin genes experienced strong positive selection and repeated gene duplications resulting in the formation of a hyperdiverse family of venom peptides. Intriguingly, the largest number of consomatins was found in worm-hunting species (>400 sequences), indicating a homologous system in annelids, another large protostome phylum. Consistent with this, comprehensive sequence mining enabled the identification of SSRP-like sequences (and their corresponding orphan receptor) in annelids and several other protostome phyla. These results established the existence of SSRP-like peptides in many major branches of bilaterians and challenge the prevailing hypothesis that deuterostome SSRPs and protostome allatostatin-C are orthologous peptide families. Finally, having a large set of predator-prey SSRP sequences available, we show that although the cone snail's signaling SSRP-like genes are under purifying selection, the venom consomatin genes experience rapid directional selection to target receptors in a changing mix of prey.

Comparative Venomics of C. flavidus and C. frigidus and Closely Related Vermivorous Cone Snails.

Cone snail venom biodiversity reflects dietary preference and predatory and defensive envenomation strategies across the ≈900 species of Conidae. To better understand the mechanisms of adaptive radiations in closely related species, we investigated the venom of two phylogenetically and spatially related species, C. flavidus and C. frigidus of the Virgiconus clade. Transcriptomic analysis revealed that the major superfamily profiles were conserved between the two species, including 68 shared conotoxin transcripts. These shared transcripts contributed 90% of the conotoxin expression in C. frigidus and only 49% in C. flavidus, which showed greater toxin diversification in the dominant O1, I2, A, O2, O3, and M superfamilies compared to C. frigidus. On the basis of morphology, two additional sub-groups closely resembling C. flavidus were also identified from One Tree Island Reef. Despite the morphological resemblance, the venom duct proteomes of these cryptic sub-groups were distinct from C. flavidus. We suggest rapid conotoxin sequence divergence may have facilitated adaptive radiation and the establishment of new species and the regulatory mechanisms facilitating species-specific venom evolution.

Venomics Reveals a Non-Compartmentalised Venom Gland in the Early Diverged Vermivorous Conus distans.

The defensive use of cone snail venom is hypothesised to have first arisen in ancestral worm-hunting snails and later repurposed in a compartmentalised venom duct to facilitate the dietary shift to molluscivory and piscivory. Consistent with its placement in a basal lineage, we demonstrate that the C. distans venom gland lacked distinct compartmentalisation. Transcriptomics revealed C. distans expressed a wide range of structural classes, with inhibitory cysteine knot (ICK)-containing peptides dominating. To better understand the evolution of the venom gland compartmentalisation, we compared C. distans to C. planorbis, the earliest diverging species from which a defence-evoked venom has been obtained, and fish-hunting C. geographus from the Gastridium subgenus that injects distinct defensive and predatory venoms. These comparisons support the hypothesis that venom gland compartmentalisation arose in worm-hunting species and enabled repurposing of venom peptides to facilitate the dietary shift from vermivory to molluscivory and piscivory in more recently diverged cone snail lineages.

Biomedical Potential of the Neglected Molluscivorous and Vermivorous Conus Species.

Within the Conidae family, the piscivorous Conus species have been a hotspot target for drug discovery. Here, we assess the relevance of Conus and their other feeding habits, and thus under distinctive evolutionary constraints, to highlight the potential of neglected molluscivorous and vermivorous species in biomedical research and pharmaceutical industry. By singling out the areas with inadequate Conus disquisition, such as the Tamil Nadu Coast and the Andaman Islands, research resources can be expanded and better protected through awareness. In this study, 728 Conus species and 190 species from three other genera (1 from Californiconus, 159 from Conasprella and 30 from Profundiconus) in the Conidae family are assessed. The phylogenetic relationships of the Conidae species are determined and their known feeding habits superimposed. The worm-hunting species appeared first, and later the mollusc- and fish-hunting species were derived independently in the Neogene period (around 23 million years ago). Interestingly, many Conus species in the warm and shallow waters become polyphagous, allowing them to hunt both fish and worms, given the opportunities. Such newly gained trait is multi originated. This is controversial, given the traditional idea that most Conus species are specialized to hunt certain prey categories. However, it shows the functional complexity and great potential of conopeptides from some worm-eating species. Pharmaceutical attempts and relevant omics data have been differentially obtained. Indeed, data from the fish-hunting species receive strong preference over the worm-hunting ones. Expectedly, conopeptides from the fish-hunting species are believed to include the most potential candidates for biomedical research. Our work revisits major findings throughout the Conus evolution and emphasizes the importance of increasing omics surveys complemented with further behavior observation studies. Hence, we claim that Conus species and their feeding habits are equally important, highlighting many places left for Conus exploration worldwide. We also discuss the Conotoxin drug discovery potentials and the urgency of protecting the bioresources of Conus species. In particular, some vermivorous species have demonstrated great potential in malaria therapy, while other conotoxins from several worm- and mollusc-eating species exhibited explicit correlation with SARS-CoV-2. Reclaiming idle data with new perspectives could also promote interdisciplinary studies in both virological and toxicological fields.

Identification of Novel Conotoxin Precursors from the Cone Snail Conus spurius by High-Throughput RNA Sequencing.

Marine gastropods of the genus Conus, comprising more than 800 species, have the characteristic of injecting worms and other prey with venom. These conopeptide toxins, highly diverse in structure and action, are highly potent and specific for their molecular targets (ion channels, receptors, and transporters of the prey's nervous system), and thus are important research tools and source for drug discovery. Next-generation sequencing technologies are speeding up the discovery of novel conopeptides in many of these species, but only limited information is available for Conus spurius, which inhabits sandy mud. To search for new precursor conopeptides, we analyzed the transcriptome of the venous ducts of C. spurius and identified 55 putative conotoxins. Seven were selected for further study and confirmed by Sanger sequencing to belong to the M-superfamily (Sr3.M01 and Sr3.M02), A-superfamily (Sr1.A01 and Sr1.A02), O-superfamily (Sr15.O01), and Con-ikot-ikot (Sr21.CII01 and Sr22.CII02). Six of these have never been reported. To our knowledge, this report is the first to use high-throughput RNA sequencing for the study of the diversity of C. spurius conotoxins.

A phylogeny-aware approach reveals unexpected venom components in divergent lineages of cone snails.

Marine gastropods of the genus Conus are renowned for their remarkable diversity and deadly venoms. While Conus venoms are increasingly well studied for their biomedical applications, we know surprisingly little about venom composition in other lineages of Conidae. We performed comprehensive venom transcriptomic profiling for Conasprella coriolisi and Pygmaeconus traillii, first time for both respective genera. We complemented reference-based transcriptome annotation by a de novo toxin prediction guided by phylogeny, which involved transcriptomic data on two additional 'divergent' cone snail lineages, Profundiconus, and Californiconus. We identified toxin clusters (SSCs) shared among all or some of the four analysed genera based on the identity of the signal region-a molecular tag present in toxins. In total, 116 and 98 putative toxins represent 29 and 28 toxin gene superfamilies in Conasprella and Pygmaeconus, respectively; about quarter of these only found by semi-manual annotation of the SSCs. Two rare gene superfamilies, originally identified from fish-hunting cone snails, were detected outside Conus rather unexpectedly, so we further investigated their distribution across Conidae radiation. We demonstrate that both these, in fact, are ubiquitous in Conidae, sometimes with extremely high expression. Our findings demonstrate how a phylogeny-aware approach circumvents methodological caveats of similarity-based transcriptome annotation.

The genome of the venomous snail Lautoconus ventricosus sheds light on the origin of conotoxin diversity.

BACKGROUND: Venoms are deadly weapons to subdue prey or deter predators that have evolved independently in many animal lineages. The genomes of venomous animals are essential to understand the evolutionary mechanisms involved in the origin and diversification of venoms. RESULTS: Here, we report the chromosome-level genome of the venomous Mediterranean cone snail, Lautoconus ventricosus (Caenogastropoda: Conidae). The total size of the assembly is 3.59 Gb; it has high contiguity (N50 = 93.53 Mb) and 86.6 Mb of the genome assembled into the 35 largest scaffolds or pseudochromosomes. On the basis of venom gland transcriptomes, we annotated 262 complete genes encoding conotoxin precursors, hormones, and other venom-related proteins. These genes were scattered in the different pseudochromosomes and located within repetitive regions. The genes encoding conotoxin precursors were normally structured into 3 exons, which did not necessarily coincide with the 3 structural domains of the corresponding proteins. Additionally, we found evidence in the L. ventricosus genome for a past whole-genome duplication event by means of conserved gene synteny with the Pomacea canaliculata genome, the only one available at the chromosome level within Caenogastropoda. The whole-genome duplication event was further confirmed by the presence of a duplicated hox gene cluster. Key genes for gastropod biology including those encoding proteins related to development, shell formation, and sex were located in the genome. CONCLUSIONS: The new high-quality L. ventricosus genome should become a reference for assembling and analyzing new gastropod genomes and will contribute to future evolutionary genomic studies among venomous animals.

Diversity of Conopeptides and Conoenzymes from the Venom Duct of the Marine Cone Snail Conus bayani as Determined from Transcriptomic and Proteomic Analyses.

Marine cone snails are predatory gastropods characterized by a well-developed venom apparatus and highly evolved hunting strategies that utilize toxins to paralyze prey and defend against predators. The venom of each species of cone snail has a large number of pharmacologically active peptides known as conopeptides or conotoxins that are usually unique in each species. Nevertheless, venoms of only very few species have been characterized so far by transcriptomic approaches. In this study, we used transcriptome sequencing technologies and mass spectrometric methods to describe the diversity of venom components expressed by a worm-hunting species, Conus bayani. A total of 82 conotoxin sequences were retrieved from transcriptomic data that contain 54 validated conotoxin sequences clustered into 21 gene superfamilies including divergent gene family, 17 sequences clustered to 6 different conotoxin classes, and 11 conotoxins classified as unassigned gene family. Seven new conotoxin sequences showed unusual cysteine patterns. We were also able to identify 19 peptide sequences using mass spectrometry that completely overlapped with the conotoxin sequences obtained from transcriptome analysis. Importantly, herein we document the presence of 16 proteins that include five post-translational modifying enzymes obtained from transcriptomic data. Our results revealed diverse and novel conopeptides of an unexplored species that could be used extensively in biomedical research due to their therapeutic potentials.

The Redox-Active Conopeptide Derived from the Venom Duct Transcriptome of Conus lividus Assists in the Oxidative Folding of Conotoxin.

The tetrapeptides Li504 and Li520, differing in the modification of the 4-trans-hydroxylation of proline, are novel conopeptides derived from the venom duct transcriptome of the marine cone snail Conus lividus. These predicted mature peptides are homologous to the active site motif of oxidoreductases that catalyze the oxidation, reduction, and rearrangement of disulfide bonds in peptides and proteins. The estimated reduction potential of the disulfide of Li504 and Li520 is within the range of disulfide reduction potentials of oxidoreductases, indicating that they may catalyze the oxidative folding of conotoxins. Conformational features of Li504 and Li520 include the trans configuration of the Cys1-Pro2/Hyp2 peptide bond with a type 1 turn that is similar to the active site motif of glutaredoxin that regulates the oxidation of cysteine thiols to disulfides. Li504- and Li520-assisted oxidative folding of α-conotoxin ImI confirms that Li520 improves the yield of the natively folded peptide by concomitantly decreasing the yield of the non-native disulfide isomer and thus acts as a miniature disulfide isomerase. The geometry of the Cys1-Hyp2 peptide bond of Li520 shifts between the trans and cis configurations in the disulfide form and thiol/thiolate form, which regulates the deprotonation of the N-terminal cysteine residue. Hydrogen bonding of the hydroxyl group of 4-trans-hydroxyproline with the interpeptide chain unit in the mixed disulfide form may play a vital role in shifting the geometry of the Cys1-Hyp2 peptide bond from cis to trans configuration. The Li520 conopeptide together with similar peptides derived from other species may constitute a new family of "redox-active" conopeptides that are integral components of the oxidative folding machinery of conotoxins.

Diversity of Conopeptides and Their Precursor Genes of Conus Litteratus.

The venom of various Conus species is composed of a rich variety of unique bioactive peptides, commonly referred to as conotoxins (conopeptides). Most conopeptides have specific receptors or ion channels as physiologically relevant targets. In this paper, high-throughput transcriptome sequencing was performed to analyze putative conotoxin transcripts from the venom duct of a vermivorous cone snail species, Conus litteratus native to the South China Sea. A total of 128 putative conotoxins were identified, most of them belonging to 22 known superfamilies, with 43 conotoxins being regarded as belonging to new superfamilies. Notably, the M superfamily was the most abundant in conotoxins among the known superfamilies. A total of 15 known cysteine frameworks were also described. The largest proportion of cysteine frameworks were VI/VII (C-C-CC-C-C), IX (C-C-C-C-C-C) and XIV (C-C-C-C). In addition, five novel cysteine patterns were also discovered. Simple sequence repeat detection results showed that di-nucleotide was the major type of repetition, and the codon usage bias results indicated that the codon usage bias of the conotoxin genes was weak, but the M, O1, O2 superfamilies differed in codon preference. Gene cloning indicated that there was no intron in conotoxins of the B1- or J superfamily, one intron with 1273-1339 bp existed in a mature region of the F superfamily, which is different from the previously reported gene structure of conotoxins from other superfamilies. This study will enhance our understanding of conotoxin diversity, and the new conotoxins discovered in this paper will provide more potential candidates for the development of pharmacological probes and marine peptide drugs.

Transcriptomic Profiling Reveals Extraordinary Diversity of Venom Peptides in Unexplored Predatory Gastropods of the Genus Clavus.

Predatory gastropods of the superfamily Conoidea number over 12,000 living species. The evolutionary success of this lineage can be explained by the ability of conoideans to produce complex venoms for hunting, defense, and competitive interactions. Whereas venoms of cone snails (family Conidae) have become increasingly well studied, the venoms of most other conoidean lineages remain largely uncharacterized. In the present study, we present the venom gland transcriptomes of two species of the genus Clavus that belong to the family Drilliidae. Venom gland transcriptomes of two specimens of Clavus canalicularis and two specimens of Clavus davidgilmouri were analyzed, leading to the identification of a total of 1,176 putative venom peptide toxins (drillipeptides). Based on the combined evidence of secretion signal sequence identity, entire precursor similarity search (BLAST), and the orthology inference, putative Clavus toxins were assigned to 158 different gene families. The majority of identified transcripts comprise signal, pro-, mature peptide, and post-regions, with a typically short (<50 amino acids) and cysteine-rich mature peptide region. Thus, drillipeptides are structurally similar to conotoxins. However, convincing homology with known groups of Conus toxins was only detected for very few toxin families. Among these are Clavus counterparts of Conus venom insulins (drillinsulins), porins (drilliporins), and highly diversified lectins (drillilectins). The short size of most drillipeptides and structural similarity to conotoxins were unexpected, given that most related conoidean gastropod families (Terebridae and Turridae) possess longer mature peptide regions. Our findings indicate that, similar to conotoxins, drillipeptides may represent a valuable resource for future pharmacological exploration.

Insertions and Deletions Play an Important Role in the Diversity of Conotoxins.

Previous studies have indicated that each conotoxin precursor has a hyperconserved signal region, a rather conserved pro region and a hypervariable mature region, and nucleotide mutations are the main driven factor. However, in this study, we made an in-depth analysis on the M-superfamily conotoxin precursors and found that the diversity of the signal, pro and mature regions are more complicated than previous findings. Different conotoxin precursors can have same signal, pro and/or mature regions, especially different conotoxin precursors with same mature region but different signal and pro regions. In addition, insertions and deletions (indels) were detected in conotoxin precursors. Indels are infrequent in the signal region but frequent in the pro and mature regions. In contrast to deletions that dominate in the pro region, insertions dominate in the mature region. The number of amino acids is crucial for the physiological functions of mature conotoxins, therefore indels, especially insertions in the mature region, play an important role in the sequence and function diversity of conotoxins.

Cone snail analogs of the pituitary hormones oxytocin/vasopressin and their carrier protein neurophysin. Proteomic and transcriptomic identification of conopressins and conophysins.

Transcriptomic analysis of cone snail venom duct tissue has permitted the identification of diverse conopressin/conophysin precursor sequences from seven distinct Conus species. Multiple precursor isoforms are present in C.monile, C.lividus and C.loroisii. Aqueous extracts of the venom duct tissue from C.monile yield a band, at ~ 15-20 kDa on SDS-PAGE. In-gel trypsin digestion, followed by mass spectrometry establishes the presence of two distinct conopressin/conophysin isoforms that differ at position 8 in the predicted conopressin nonapeptide sequence. Mass spectrometric analysis of aqueous extracts revealed the presence of four conopressin related peptides, whose sequences could be deduced from MS/MS fragmentation patterns. The four sequences determined in this study are CFIRNCPKG*, CFIRNCPEG*, CFIRNCPK* and CFIRNCPE* (∗ indicates amide), which were further confirmed by comparison with chemically synthesized peptides. A conophysin with a mass of 9419.7 Da was also detected, corresponding to one of the isoforms revealed by the transcriptome data. Complete conservation of fourteen Cys residues and the key residues involved in peptide hormone binding is established by comparison of conophysin sequences, with the crystallographically characterized sequence of bovine neurophysin, in complex with vasopressin. A survey of available sequences for oxytocin/vasopressin peptides in both vertebrates and invertebrates establishes the conopressins as a distinct group in this family. C-terminal amidated, truncated conopressin analogs may arise by alternate post-translational processing.